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Abstract 

We present the full 0(a) electroweak radiative corrections to single Higgs production in e + e~ . 
This takes into account the full one-loop corrections as well as the effects of hard photon radia- 
tion. We include both the fusion and Higgs-strahlung processes. The computation is performed 
with the help of GRACE-loop where we have implemented a generalised non-linear gauge fixing 
condition. The latter includes 5 gauge parameters that can be used for checks on our results. 
Besides the UV, IR finiteness and gauge parameter independence checks it proves also powerful 
to test our implementation of the 5-point function. We find that for a 500GeV machine and a 
light Higgs of mass 150GeV, the total 0(a) correction is small when the results are expressed 
in terms of «qed- The total correction decreases slightly for higher energies. For moderate 
centre of mass energies the total 0(a) decreases as the Higgs mass increases, reaching —10% for 
Mh = 350GeV and y/s = 500GeV. In order to quantify the genuine weak corrections we have 
subtracted the universal virtual and bremsstrahlung correction from the full 0(a). We find, for 
Mh = 150GeV, a weak correction slowly decreasing from —2% to —4% as the energy increases 
from = 300GeV to = ITeV after expressing the tree- level results in terms of G^. 
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1 Introduction 



Uncovering the mechanism of symmetry breaking is one of the major tasks of the high 
energy colliders. Most prominent is the search for the Higgs particle. Although the 
LHC should not miss this particle even if it weighed up to ITeV, precision measure- 
ments on the Higgs properties will only be conducted in an e + e~ collider. There are 
two important mechanisms for Higgs production in e + e _ . The Higgs-strahlung process, 
e + e~ — ► ZH and the PU-fusion process, e + e~ — » v e v e H . The former is the dominant 
one at small (in the LEP2 range say) to moderate energies but decreases rather fast with 
energy. At TeV energies the H^-fusion process dominates by far for Higgs masses up to 
ITeV. Even at 500GeV this t-channel process is dominant for Higgs masses in the range 
preferred by the indirect electroweak precision data|I] and remains an important com- 
ponent for all Higgs masses at energies of the linear collider. Tree-level computations of 
Higgs production are rather well under control]^, [J, including interference of the W fusion 
process with the Higgs-strahlung process. Note however that the complexity of the process 
e + e _ — > vvH precludes a full analytic result for the total cross section even at tree-level, 
although the differential cross section can be cast in a very compact formf^, [|. Full 
radiative corrections for Higgs-strahlung have been considered by a number of groupsQ, 
while a proper one-loop treatment of the fusion process is still lacking despite the im- 
portance of the process for the linear collider physics program. Some recipes have been 
suggested § to include parts of the radiative corrections to the fusion process but con- 
sidering the domain of validity of these approximations y/s, Mh <C 2m t , they are expected 
not to be precise for the interesting range of Higgs masses (preferred by the latest precision 
measurements ||) and next collider energies. One-loop contribution to the HWW vertex 
has been considered on the basis that it might constitute a good approximation for the 
fusion process^, but it rests to see how well this approximation fares in comparison of 
the full calculation. Very recently one-loop radiative corrections to this process have been 
investigated within the minimal supersymmetric model but again by only taking into ac- 
count the contribution of the fermions and sfermions to the H/hWW vertex]?], |J. It is 
the aim of this letter to summarise the results of the full radiative corrections to single 
Higgs production in e + e~ , including both the fusion and Higgs-strahlung processes in the 
SM. (Standard Model). We include both the virtual and soft corrections as well as the 
hard photon radiation. A longer paper will detail our computation and results and will 
look into the issue of finding approximations to the full result[[ 

1 Preliminary results have been presented at the Workshop RADCOR2002§. At this meeting the 
FIRCLAjlO group exposed their plans and techniques, different from ours concerning Feynman integration, 
for tackling the calculation of this process. While finalising this letter we also learnt of a calculation by 
A. Denner, S. Dittmaier, M. Roth and M. Weber, in preparation. 
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A standard hand calculation using the usual techniques could hardly be attempted for 
such 2^3 processes at one-loop. Considering the ever increasing power of computers, the 
possibility of parallelisation and the fact that the whole procedure of perturbation theory 
consists of algorithms that can be directly translated on a computer it seems that most, if 
not all, complex calculations in high-energy physics can be automated. This is especially 
true for electroweak processes where various scales and masses enter the calculations. 



GRACE-loop []TT[ from which our results are derived is such a program. GRACE JT2|], the 
tree-level component of the system, has been tested and heavily used for tree-level cross 



sections up to 6-fermions in the final state[13|. GRACE-loop has been exploited and 
checked thoroughly for a variety of 2 —>■ 2 processes in the electroweak theory [14|. The 
system which requires as input, a model file that describes all the interaction vertices 
derived from a particular Lagrangian can generate all the necessary Feynman graphs 
together with their codes so that matrix elements can be generated before being processed 
for the calculation of the cross section and event generation. For loop processes, there is 



a symbolic manipulation stage (either F0RM[15] or REDUCE p"6|l ) that handles all the Dirac 
and tensor algebra in n-dimension for all the interference terms between tree-level and 
1-loop diagrams and automatically applies the Feynman trick for the propagator. This 
is then passed to a module that contains two libraries for the loop integration containing 
the FF package [jT7H as well as an in-house numerical code. The system together with the 
one-loop renormalisation program is described in detail in As far as the calculation 
of one-loop processes is concerned, a series of powerful tests are implemented in the code 



as described in 14 and as will be presented below for e + e — > uvH . 



2 Tree-level results, setting-up the loop calculation 

Our input parameters for the calculation of e + e~ — > vvH are the following. Throughout 

we expressed our results in terms of the fine structure constant in the Thomson limit 

or x = 137.0359895 and the Z mass Mz = 91.1876GeV. Our on-shell renormalisation 

program uses My/ as input parameter, nonetheless our numerical value of Mjy is derived 

through Ar]l8|]0. M\y thus changes as a function of Mu- For the the lepton masses we 

take m e = 0.510999MeV, m M = 105.6584MeV and m T = 1.777GeV. For the quark masses 

beside the top mass M t = 174GeV, we take the set M u = M d = 58MeV, M s = 92MeV, 

M c = 1.5GeV and M b = 4.7GeV. With these values we calculate Aa(M z ) = 0.059258. 

With this we find for example that M w = 80.3767GeV for M H = 150GeV and M w = 

80.3158GeV for Mu = 350GeV. Especially for the Higgs-strahlung subprocess we require 

2 We include NLO QCD corrections and two- loop Higgs effects. We take a s (Af§) = .118 together with 
G At = 1.16639 x l(T 5 GeV~ 2 



2 



a Z-width. We have taken a constant fixed Z-width, Y z = 2.4952GeV. Unless when 
otherwise stated our results refer to the full e + e~ — > vvH , summing over all three types 
of neutrinos with, for electron neutrinos, the effect of interference between fusion and 
Higgs-strahlung. 

We have checked that our tree-level results are in very good agreement with those in 
after expressing them in terms of G^. However since we are considering the effect of 
radiative corrections, within our scheme we prefer showing all our results using a. We will 
only comment on the scheme at the end of this letter. We find for example for y/s = 
500GeV and Mh = 350GeV that cr treei totai = 4.606fb at tree- level for the contribution of 
all three neutrinos. In an attempt to separate the different contributions to single Higgs 
production, we will refer to the s-channel as given by a s = 3 x cr(e + e~ — > PyV^H). The 
bulk of this contribution is given by cr(e + e~ — > ZH) x Bz-^ mv . , Bz->mv. — 20%. We will 
define the t-channel as at = Ctotai — a s , this implicitly means that the interference term is 
included in this contribution. These definitions will be carried over to the one-loop case 
as well. 

In Fig. H we have also included the tree-level cross section. They are shown as a 
function of the centre-of-mass energy for a light Higgs of mass Mu = 150GeV as well as a 
function of the Higgs mass at a centre-of-mass of 500GeV. All integration over phase space 



are done with the help of BASES, see|12|. These figures clearly show the importance of the 
t-channel contribution pointed out in the introduction. For a low Higgs mass of 150GeV, 
although the s-channel still dominates at y/s = 300GeV, very quickly at y/s = 500GeV 
it is the t-channel that dominates. For the latter energy as the Higgs mass increases, 
the t-channel contribution drops much quickly than the s-channel, where both merge 
around M H = 380GeV to M H = 390GeV, but very quickly around the ZH threshold 
for Mh ~ 408GeV, the s-channel drops precipitously leaving the t-channel as the sole 
contribution for the whole process. 

Neglecting all Goldstone-electron coupling (proportional to the electron mass), one 
has at one-loop 249 diagrams (for i> e u e H, and 146 for v^v^H) including 15 pentagons 
(5-point functions) compared to only 2 diagrams at tree-level, one for the fusion process 
and one for the Higgs-strahlung. Keeping the electron Yukawa coupling one has a total of 
1350 diagrams (98 pentagons corresponding to 5-point functions) at one-loop. Although 
in running our program to derive cross sections we only use the set of 249 diagrams, to 
perform our extensive checks especially those of gauge-parameter independence, at the 
level of the differential cross section, we keep the full set of 1350 diagrams. It is impossible 
to show all the contributing diagrams here. They may be downloaded or visualised at 
this location [19|]. All these diagrams are generated and drawn by gracef ig the Feynman 
diagrams generator of GRACE . A representative selection of diagrams is shown in Fig. [L[ 
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Graph 87 Graph 249 Graph 486 




produced by GRACEFIG 



Figure 1: A small selection of different classes of loop diagrams contributing to e + e~ — > 
vvH . We keep the same graph numbering as that produced by the system. Graph 1312 
belongs to the corrections from self- energies, here both the virtual and counterterm contri- 
butions are generated and counted as one diagram. Graph 87 shows a vertex correction. 
Both graphs can be considered as resonant Higgs-strahlung contributions. Graph 249 rep- 
resents a box correction, it is a non resonant contribution but applies also to the v T 
channels. Graph 486 is also a box correction which is non resonant and applies only to u e . 
Graph 541 and Graph 565 are typical bosonic and fermionic corrections to the WWH 
vertex for the fusion process. Graph 846 shows a pentagon correction that also applies to 
fj, and r neutrinos, this again can be considered as a non-resonant contribution. Graphs 
827 and 828 are pentagons that only contribute to e + e~ — > v e v e H . 
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The results of the calculation are checked by performing three kinds of tests at some 
random points in phase space. For these tests to be passed one works in quadruple 
precision. We first check the ultraviolet finiteness of the results. This test applies to the 
whole set of the virtual one-loop diagrams. In order to conduct this test we regularise any 
infrared divergence by giving the photon a fictitious mass (we set this at A = 10~ 15 GeV). 
In the intermediate step of the symbolic calculation dealing with loop integrals (in n- 
dimension), we extract the regulator constant Cjjv — V £ — 1e + log47r, n = 4 — 2e and 
treat this as a parameter. The ultraviolet finiteness test gives a result that is stable over 
30 digits when one varies the dimensional regularisation parameter Cuv- This parameter 
could then be set to in further computation. The test on the infrared finiteness is 
performed by including both loop and bremsstrahlung contributions and checking that 
there is no dependence on the fictitious photon mass A. We find results that are stable 
over 23 digits when varying A. An additional test concerns the bremsstrahlung part. It 
relates to the independence in the parameter k c which is a soft photon cut parameter 
that separates soft photon radiation and the hard photon performed by the Monte-Carlo 
integration. 

Gauge parameter independence of the result is performed through a set of five gauge 
fixing parameters. For the latter a generalised non- linear gauge fixing condition f2(| has 
been chosen. 

C GF = -^H - ie & A n - igcwPZ^)W» + + £ w 9 -(v + ~5H + iR X z)X + \ 2 
Kw 1 

1 -{d.Z + £z-?-(v + £H) X3 ) 2 - ^-{d.Af . (2.1) 



2£ z v ^2c w y 26 

The x represent the Goldstone. We take the 't Hooft-Feynman gauge with £w = Cz = 
£a = 1 so that no "longitudinal" term in the gauge propagators contributes. Not only this 
makes the expressions much simpler and avoids unnecessary large cancelations, but it also 
avoids the need for high tensor structures in the loop integrals. The use of five parameters 
is not redundant as often these parameters check complementary sets of diagrams. For 
example the parameter (3 is involved in all diagrams containing the gauge WW Z and their 
Goldstone counterpart, whereas a checks WWj and 5 is implicitly present in WWH . For 
each parameter of the set ( = (a, j3,5, k, e) the first check is made while freezing all other 
four parameters to 0. We have also made checks with two parameters non-zero. In 
principle checking for 2 or 3 values of the gauge parameter should be convincing enough. 
We in fact go one step further and perform a complete gauge parameter independence. 
To achieve this we generate for each non-linear gauge parameter (, the values of the 
loop correction to the total differential cross section as well as the contribution of each 
one-loop diagram contribution for the five values ( = 0,±1,±2. The one-loop diagram 
contribution from each loop graph g, is defined as 
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da g = da g (() = & {jl°° v ■ T tree f ) . (2.2) 

T tree is the tree-level amplitude summed over all tree-diagrams|. T g loop is the one-loop 
amplitude contribution of the one-loop diagram g. A rapid look at the structure of the 
Feynman rules of the non-linear gauge leads one to conclude that for e + e~ — > vvH each 
contribution is a polynomial of (at most) third degree in the gauge parameter and thus, 
that each contribution, da g may be written as 



d* g = daf + Cda« + C 2 daf + (Ma^ , (2 . 3 ) 

For each contribution da g , it is a straightforward matter, given the values of da g for 
the five input £ = 0, ±1, ±2, to reconstruct dcr^ 0,1 ' 2,3 -'. This is what we do. In fact for each 
set of parameters we automatically pick up all those diagrams that involve a dependence 
on the gauge parameter. The number of diagrams in this set depends on the parameter 
chosen. In some cases a very large number of diagrams is involved. For the process at 
hand this occurs with the parameter 5 where about 500 diagrams are involved in the 
check. 

We then verify that the differential cross section is independent of ( 



da = £da 9 = £da(°), (2.4) 



and therefore that 



^9 UL ff 

Max 9 (|d4' ; |) 

vanishes. 



sum, = - 9 _\ t) , i = 1,2,3, (2.5) 



As seen from Table [TJ agreement within 20 to 30 digits is observed. This agreement 
gets better if one gives the electron mass a higher value, say lGeV. The gauge parameter 



3 Therefore the tree-level amplitude does not depend on any gauge parameter. For the process at hand 
nonetheless, some individual tree diagrams depend on the gauge parameter 8 and e giving extremely small 
contributions proportional to the electron mass. Our numerical procedure to isolate the gauge parameter 
dependence detects these tiny variations, replacing T tree < by any tree-level diagram instead of the sum, 
one is able to differentiate between a variation in da g due to the loop diagram or the residual one from 

rj-trt-e f 
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# graphs 


sum 3 


sum 2 


sumi 


a 


149 




1Q -28 


1Q -30 




314 




10- 31 


1Q -23 


5 


477 (1059) 


10 -20 


10 -20 


1Q -26 


K 


122 




1Q -23 


1Q -23 


6 


128 (132) 




icr 21 


1Q -30 



Table 1: Numerical size of surrii for each non-linear gauge parameter. '# of graphs' means 
the number of graphs that contributes to each sum depending on the gauge parameter, the 
number between bracket include the residual dependence from T tree * . — means that no 
diagram is involved. 



dependence check not only tests the various components of the input file (correct Feynman 
diagrams for example) but also the symbolic manipulation part and most important of all 
the correctness of all the reduction formulae and the proper implementation of all the Ap- 
point functions. This is quite useful when one deals with 5-point functions as is the case at 
hand. Talking of parametric integrals all tensor reductions are done following the standard 
procedure and then passing the scalar integrals to the FF package fi~7jl or to our own 



specially optimised routines when photon exchange is involved. The pentagon integrals 



are expressed in terms of boxes as is now standard |21j, our procedure is outlined in the 
Appendix. We work in the on-shell renormalisation scheme closely following ||22|| . Apart 
from masses and couplings, renormalisation is also carried for the fields. In particular 
we also require the residues of the renormalised propagators of all physical particles to 
be unity. As known|4[], this procedure leads to a (very sharp) threshold singularity in 
the wave function of the Higgs at the thresholds corresponding to M# = 2Mw,2Mz- 



Solutions to smooth this behaviour |[23||, like the inclusion of the finite width of the W and 
Z, do exist but we have not implemented them yet in the present version of GRACE- loop . 
Therefore when scanning over M# it is sufficient to avoid these regions within lGeV 
around the thresholds. 

As a separate check on our implementation we have also computed H — > WW 
and e + e~ — ► ZH, after tuning our parameters we find excellent agreement with the 
literature [§J |J. 
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3 Results 



3.1 Full O(o) results 

The results we show here include all 3 neutrino species. We first discuss the full 0(a) 
which includes the hard bremsstrahlung part. The separation between the s-channel and 
the t-channel is done in the same way as with the separation done at tree-level. It should 
be noted that in the one-loop diagrams that contribute to e + e~ — > v^v^H, and which 
we classify as s-channel, there are diagrams which can not be deduced from the one-loop 
corrections to the s-channel 2-^2 process e + e~ — > ZH. Graph 249 in Fig.|l] is one 
example. The relative correction for all three contributions is defined as 5o(a) = — 1 
and will be referred to as the full one-loop 0(a) correction. 

One important remark is that the overall correction in the s-channel and i-channel 
are quite different. For a light Higgs of mass 150GeV the correction to the s-channel 
Higgs-strahlung contribution is positive for practically all centre of mass energies of the 
next linear collider, see Fig. ^. It rises rather sharply to reach about +20% for a centre-of 
mass energy of ITeV, however as will be argued below, the bulk of these large corrections 
are due to virtual and real QED corrections. Moreover in regions where these corrections 
are large, the Higgs-strahlung contribution is rather small. On the other hand the total 
correction in the t-channel for a small Higgs mass of 150GeV is negative throughout the 
range y/s = 350GeV to ITeV, and is almost constant past ^/s = 500GeV, reaching about 
—2%. Combining these two contributions, we see that the full correction to the whole 
process also remains small for a small Higgs mass. In fact at 500GeV it is almost at 
its lowest of about 0.2%. This is an accidental cancellation between the contributions 
of both the t-channel and the s-channel at this energy. We have also studied the Higgs 
mass dependence of the corrections. First of all note that all our results capture the sharp 
spikes for Mh = 2Myy,2Mz, the top threshold is also visible when we plot the relative 
corrections. For 500GeV the s-channel 0(a) correction decreases for higher Higgs masses, 
eventually turning negative with a value —6% for Mh = 350GeV. It then drops rather 
sharply to reach as much as —30% at the ZH threshold, past which this cross section is 
completely negligible. This behaviour is largely due to QED corrections and is driven by 
the kinematics of the two-body e + e~ — ► ZH . At 500GeV the correction in the i-channel 
contribution remains negative for all Higgs masses that we considered, i.e. in the range 
115 - 450GeV. It drops steadily from about -2% for M H = 115GeV to about -10% at 
350GeV close to the top pair threshold. It then increases up to the ZH threshold before 
dropping sharply around the ZH threshold. Most of the large corrections are due to QED 
corrections. 
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150 2 2 5 3 3 5 4 4 50 150 200 250 500 550 400 450 



M H M H (GeV) 

Figure 2: The two figures in the first row show the cross section as a function of centre-of- 
mass energy for a light Higgs of mass Mh = 150 GeV. We show the s-channel, t-channel 
and the sum of these (total) cross sections as defined in the text. Both the tree-level 
(dashed lines) and the full one-loop correction (full lines) are shown. In the second panel 
we show the relative correction in per-cent. In the second row, the dependence of the cross 
section as a function of the Higgs mass at a centre- of-mass of 500GeV is shown. 
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3.2 Extraction the QED corrections 



As is known large QED corrections require a higher order treatment. In order to quantify 
the effects of the genuine weak corrections, one could try subtract these QED corrections. 
This can be done rather easily for the s-channel contribution, where the correction can 
be readily extracted from the electromagnetic correction to the eeZ vertex and the soft- 
photon bremsstrahlung part. Indeed our computation produces at an intermediate stage 
the result including the soft bremsstrahlung correction, that is before the inclusion of hard 
photons. The cut on the photon energy, k c , has been taken sufficiently small, k c = O.lGeV. 
These corrections without hard bremsstrahlung include thus the QED virtual and soft 
bremsstrahlung (which depend on k c ) as well as the genuine weak correction to the process. 
For this s-channel process the latter QED corrections are given by the universal soft 
photon factor that leads to a relative correction 



= v ( (Le - 1} ln k + l Le + t ~ x ) ' Le = Hs,ml) ■ (3 - 6) 

where m e is the electron mass and Eb the beam energy s = ^E\. Subtracting this contri- 
bution from our k c dependent (numerical) result reproduces the genuine weak correction, 

^W, s — channel • 

To quantify in an unambiguous way the effect of the weak correction in the t-channel 
we have also subtracted this universal factor from the full 0(a) correction. It can 

be shown that the leading (infrared and collinear) contributions are given by the universal 
factor, Eq. p.6|HT4| . This procedure also paves the way to a resummation of these large 
QED factors for the full process which is conducive to a Monte-Carlo implementation 
as could be done for instance through a QED parton shower |24|j . This will be treated 
elsewhere. Coming back to the weak corrections, we will denote by 5\y,t-channei and 5w,totai 
the weak correction for the t-channel and the full process based on the subtraction of the 
universal QED factor in Eq. |3.6|. 



For a light Higgs mass (M# = 150GeV), the weak correction for the s and t-channels 
have a different behaviour as the energy increases, see Fig. [| The former varies from 
about 6% at <Js = 300GeV to -2.5% at y/s = ITeV. Past 400GeV where it dominates, 
the weak correction to the t-channel varies rather slowly from 7% at y/s = 400GeV to 
about 5% at ^/s = ITeV. The dependence of the weak corrections on the Higgs mass for 
a moderate centre-of-mass energy, y/s = 500GeV, reveals that up to the ZH threshold 
these corrections increase with the Higgs mass (most probably due to Mjj terms from 
the Higgs self-coupling as in e + e~ — > ZH), apart from the clearly visible spikes at the 
W, Z and the top thresholds. Apart from the drop in the t-channel contribution around 
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the ZH threshold, the weak correction in the t-channel picks up again and as expected 
merges with the correction to the full process. 




Figure 3: Relative weak corrections as defined in the text, for the t-channel (b~w,t- channel ) , 
s-channel (b~w,s-channei)and the whole process (b~w,totai) ■ We also show the full 0(a) cor- 
rection for the whole process in the first panel. Also shown is the weak correction for the 
full process expressed in the scheme (S^ total ), see text. 



3.3 Expressing the weak corrections in terms of 

Expressing the corrections in the G a scheme or in other words had we expressed our tree- 
level results in terms of G n , thus subtracting some universal weak corrections (essentially 
fermionic contributions) affecting two-point functions, we can have a quantitative measure 
of the non-universal weak radiative correction specific to this process. We thus define 
for the s-channel and t-channel contributions, these weak corrections as 5^ = 5w — 
3Ar. Let us briefly summarise our findings for M# = 150GeV where with our input Ar 
contributes about 3% (the leading Higgs mass dependence in Ar is logarithmic). For the 
full contribution with all three neutrinos we find 5^ to be slowly varying (with exactly 
the same "slope" as 5w,totai in Fig. 0), from about —2% to about —4% in the energy range 
from y/s = 300GeV to \fs = ITeV. These genuine weak corrections remain therefore well 
contained in the full process, but in view of the precision of the e + e~ machine they must 
be taken into account. Applied to the s-channel with Mn = 150GeV, the corrections 
with G u as an input, are moderate for energies up to 400GeV but they quickly decrease 
below about —12% at ITeV. Such behaviour had been observed in e + e~ — > ZH^. 
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This is another manifestation of the failure of the scheme to properly describe the 
weak corrections for such processes at high energies. For example it is known that in 
e + e~ — > ZH the contribution of boxes is important. We do not attempt in this letter 
to make a thorough investigation of the different loop contributions to e + e~ — > uuH, for 
example the fermionic and bosonic contributions. We leave this to a further study. This 
could be interesting in order to devise reliable approximations based on a small subset 
compared to the large number of contributions for such a complex process. For example, 
very recently, the fermionic contributions and especially the effect of the third generation 
have been investigated in || and @ with differing results. It could be interesting to see 
how well these contributions can reproduce the full result. We also do not report here 
on how the distributions in the Higgs variables are affected by the radiative corrections. 
We have briefly discussed this in a previous note|| and leave the full discussion for a 
forthcoming paper. 



4 Conclusions 

We have calculated with the help of GRACE-loop the full radiative corrections including 
hard photon radiation to the important Higgs discovery channel at a future high energy 
e + e~ machine, e + e~ — > vvH . Apart from the usual checks on the ultraviolet and infra- 
red fmiteness of the result, we have performed tests on the gauge parameter independence 
of the results. To this end we have relied on a generalised non-linear gauge fixing con- 
dition where one has control over five independent gauge parameters. For a light Higgs 
of mass 150GeV for energies ranging from 300GeV to ITeV we find a modest total 0(a) 
correction which is within ±2%, being negligible at 500GeV (2 per-mil). We have also 
studied the Higgs mass dependence at ^/s = 500GeV. For example with M H = 350GeV 
we find a larger 0(a) negative correction of about —10%. In order to quantify the weak 
correction we have subtracted the universal QED virtual and bremsstrahlung corrections. 
In the energy range = 300GeV to y/s = ITeV we find, for M# = 150GeV, that for 
the full process the correction ranges from +7% to +5% when the tree-level is expressed 
in terms of a. Further investigations and details on this important process are left to a 
forthcoming publication. 
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A Appendix 



Five point functions are calculated as linear combinations of four point functions [21]. Our 
method is based on an identity suitable for the Feynman parameter integration, which is 
similar to the one described in[|25| . 



A five point function is expressed as 

dH N{1) 



(A. 1) 



(2tt) 4 AAAAA' 

where I is the loop momentum and N(l) is a polynomial of I 2 and inner products of I with 
other four-vectors. The denominators of propagators are defined as 

D = I 2 -m 2 = l 2 + X , 

A = {l + rtf-m 2 = l 2 + 2l.r i + X l , z = l,...,4. (A. 2) 

We take a set r, (i = 1,...,4) of linearly independent momenta. The latter form a 
basis for vectors in 4-dimensional space. Therefore with the Gram matrix Aij = Ti.Tj one 
has the following identity 

r = £ ,f,l,;ry =► l 2 = Y: /•'•,!,//./> (A. 3) 

i,j=l i,j=l 



Combining this identity with Eq. flA. 2 ) we obtain 



4 4 

i = EK+E' T iM fl <" ( A - 4 ) 



where 



a=0 i=l 



a i = IEA?(^-Xo), (A. 5) 
4 4 

a o = -x~Y. a i' ( A - 6 ) 

^ i=l 

Kk = ~A7 k \ (A. 7) 

4 



bo,k = - E Kk, ( A - 8) 

A = 4X + ]T(X J -X )4- 1 (X J -Xo). (A. 9) 
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This immediately shows that the five point tensor integral can be reduced to 5 box inte- 
grals. 

Now we introduce the Feynman parameters. It is easy to see that 

h = I W\i / A - E E U + E * W) *(*«)■ (A. 10) 

V / 17 A=0 0=0 ^ " V i=l / 

Making a shift in the loop momentum , I —> I — t t = X) i=1 %i r i, so as to eliminate 
linear terms in the loop momentum in D, we obtain our reduction formula 



15 = / r^\i I II <fo A <?(i - E x p) N<yl m - E f a « +E Lr i fe «,i > ) <y(sa) 

\ / J A=0 /3=0 ^ a=0 V i=l / 

4 4 

a = a - 2/ A , a, = a, , D = I 2 + ^ X a x a - ^ Ay x&j. (A. 11) 

a = M=l 

For the scalar pentagon, N(l) = 1, only a a in the previous equation contributes. 
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